Implementing a Partitioned 2-page Book Embedding 
Testing Algorithm* 



Patrizio Angelini 1 , Marco Di Bartolomeo 1 ' 2 , and Giuseppe Di Battista 

1 Dip. di Informatica e Automazione, Roma Tre University, Italy 

2 Italian Inter-University Computing Consortium CASPUR, Italy 

{angelini, gdb}@dia . uniroma3 . it , m . dibartolomeoScaspur . it 



Abstract. In a book embedding the vertices of a graph are placed on the "spine" 
of a "book" and the edges are assigned to "pages" so that edges on the same 
page do not cross. In the PARTITIONED 2-PAGE BOOK EMBEDDING problem 
egdes are partitioned into two sets E\ and E2, the pages are two, the edges of Ei 
are assigned to page 1, and the edges of Ei are assigned to page 2. The problem 
consists of checking if an ordering of the vertices exists along the spine so that the 
edges of each page do not cross. Hong and Nagamochi [13] give an interesting 
and complex linear time algorithm for tackling PARTITIONED 2-PAGE BOOK 
EMBEDDING based on SPQR-trees. We show an efficient implementation of this 
algorithm and show its effectiveness by performing a number of experimental 
tests. Because of the relationships [13] between PARTITIONED 2-PAGE BOOK 
EMBEDDING and clustered planarity we yield as a side effect an implementation 
of a clustered planarity testing where the graph has exactly two clusters. 



1 Introduction 

In a book embedding [14] of a graph the vertices are placed on the "spine" of 
a "book" and the edges are assigned to "pages" so that edges on the same page 
do not cross. A rich body of literature witnesses the interest of the scientific 
community for book embeddings. See, e.g., [3, 16]. 

Several constrained variations of book embeddings have been studied. In [15] 
the problem is tackled when in each page the number of edges incident to a ver- 
tex is bounded. In [10] the graph is directed upward planar and the order of 
the vertices on the spine must be consistent with the orientation of the edges. 
Hong and Nagamochi [13] provide a linear time algorithm for a problem called 
PARTITIONED 2-PAGE BOOK Embedding (P2BE). In the P2BE problem the 
egdes of an input graph G(V, E\, E2) are partitioned into two sets E\ and E2, 
the pages are just two, the edges of E\ are assigned to page 1, and the edges of 
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Ei are assigned to page 2. The problem consists of checking if an ordering of 
the vertices exists along the spine so that the edges of each page do not cross. 

In [13] the P2BE problem is characterized in terms of the existence of an 
embedding of G allowing to build a variation of the dual graph containing a 
particular Eulerian tour. The existence of such an embedding is tested exploiting 
SPQR-trees [8] for biconnected components and BC-trees for connected ones. 

In this paper we discuss an implementation of the algorithm in [13]. To 
efficiently implement the algorithm we faced the following problems: (i) One 
of the key steps of the algorithm requires the enumeration and the analysis of 
all the permutations of a set of objects. Even if the cardinality of the set is 
bounded by a constant this may lead to very long execution times. We restated 
that step of the algorithm avoiding such enumerations, (ii) Some steps of the 
algorithm are described in [13] at a high abstraction level. We found how to 
efficiently implement all of them, (iii) The algorithm builds several embeddings 
that are tested for the required properties only at the end of the computation. Our 
implementation considers only one embedding that is greedily built to have the 
properties. We performed experiments over a large set of suitably randomized 
graphs. The experiments show quite reasonable linear execution times. 

The algorithm in [13] is interesting in itself, since book embedding problems 
are ubiquitous in Graph Drawing. However, it is even more appealing because it 
yields [13] almost immediately a linear time algorithm for the following special 
case of clustered planarity testing. A planar graph G(Vi, V2, E) whose vertices 
are partitioned into two sets (clusters) V\ and V2 is given. Is it possible to find 
a planar drawing for G such that: (i) each of V\ and V2 is drawn inside a sim- 
ple region, (ii) the two regions are disjoint, and (iii) each edge of E crosses 
the boundary of a region at most once? Using the terminology of Clustered 
Planarity, this is a clustered planarity testing for a flat clustered graphs with ex- 
actly two clusters. References on clustered planarity can be found, e.g., in [9, 
5]. Hence, we yield, as a side effect, an implementation of such special case 
of clustered planarity testing. An alternative algorithm for the same clustered 
planarity problem has been proposed in [2, 1]. 

The paper is organized as follows. In Section 2 we give preliminaries. In 
Section 3 we outline the algorithm. Section 4 discusses how to search an em- 
bedding with the desired features and Section 5 gives further implementation 
details on the search. In Section 6 we describe our experiments. Section 7 gives 
concluding remarks. 

2 Preliminaries 

In this section we give preliminary definitions that will be used in the paper. 



2.1 Planarity 



A planar drawing of a graph is a mapping of each vertex to a distinct point of 
the plane and of each edge to a simple Jordan curve connecting its endpoints 
such that the curves representing the edges do not cross but, possibly, at com- 
mon endpoints. A graph is planar if it admits a planar drawing. Two drawings 
of a graph are equivalent if they determine the same circular ordering around 
each vertex. An embedding is an equivalence class of drawings. A planar draw- 
ing partitions the plane into topologically connected regions, called faces. The 
unbounded face is the outer face. 

2.2 Connectivity and SPQR-trees 

A graph is connected if every two vertices are joined by a path. A graph G is 
biconnected (triconnected) if removing any vertex (any two vertices) leaves G 
connected. 

To handle the decomposition of a biconnected graph into its triconnected 
components, we use SPQR-trees (see [7, 8, 12]). 

A graph is st-biconnectible if adding edge (s, t) to it yields a biconnected 
graph. Let G be an st-biconnectible graph. A separation pair of G is a pair 
of vertices whose removal disconnects the graph. A split pair of G is either a 
separation pair or a pair of adjacent vertices. A maximal split component of G 
with respect to a split pair {u, v} (or, simply, a maximal split component of 
{u, v}) is either an edge (u, v) or a maximal subgraph G' of G such that G' 
contains u and v, and {u, v} is not a split pair of G'. A vertex w ^ u,v belongs 
to exactly one maximal split component of {u, v}. We call split component of 
{u, v} the union of any number of maximal split components of {u, v}. 

We assume consider SPQR-trees that are rooted at one edge of the graph, 
called the reference edge. 

The rooted SPQR-tree T of a biconnected graph G, with respect to a refer- 
ence edge e, describes a recursive decomposition of G induced by its split pairs. 
The nodes of T are of four types: S, P, Q, and R. Their connections are called 
arcs, in order to distinguish them from the edges of G. 

Each node \i of T has an associated st-biconnectible multigraph, called the 
skeleton of fi and denoted by skel(/i). Skeleton skel(/i) shows how the children 
of \i, represented by "virtual edges", are arranged into \i. The virtual edge in 
skel(/i) associated with a child node v, is called the virtual edge of v in skel(ix). 

For each virtual edge e, of skel(//), recursively replace e, with the skeleton 
skel(/ij) of its corresponding child fa. The subgraph of G that is obtained in this 
way is the pertinent graph of \i and is denoted by pert( / u). 



Given a biconnected graph G and a reference edge e = (u',v'), tree T is 
recursively defined as follows. At each step, a split component G*, a pair of 
vertices {u, v}, and a node v in T are given. A node /i corresponding to G* 
is introduced in T and attached to its parent v. Vertices u and v are the poles 
of /i and denoted by u{fj) and v(p), respectively. The decomposition possibly 
recurs on some split components of G* . At the beginning of the decomposition 
G* = G — {e}, {u, v} = {u', v'}, and v is a Q-node corresponding to e. 

Base Case: If G* consists of exactly one edge between u and v , then /i is a 
Q-node whose skeleton is G* itself. 

Parallel Case: If G* is composed of at least two maximal split components 
Gx, . . . , Gfe (k > 2) of G with respect to {u, v}, then /i is a P-node. Graph 
skel(/i) consists of k parallel virtual edges between u and v, denoted by 
e±, . . . , e/j and corresponding to G\, . . . , Gk, respectively. The decomposi- 
tion recurs on G\, . . . , Gk, with {u, v} as pair of vertices for every graph, 
and with /i as parent node. 

Series Case: If G* is composed of exactly one maximal split component of 
G with respect to {u, v} and if G* has cutvertices ci, . . . ,Cfc_i (k > 2), 
appearing in this order on a path from u to v, then fi is an S-node. Graph 
skel(/i) is the path e±, . . . , e^, where virtual edge connects Cj_i with q 
(i = 2, . . . , k — 1), ei connects u with ci, and e& connects Cfc_i with 
The decomposition recurs on the split components corresponding to each of 
ei, e2, ■ ■ ■ , efc_i, efc with /i as parent node, and with {u, c\}, {c\, C2}, . . . , 
{ c fe-2) Cfc_i}, {ck_i,v} as pair of vertices, respectively. 

Rigid Case: If none of the above cases applies, the purpose of the decomposi- 
tion step is that of partitioning G* into the minimum number of split compo- 
nents and recurring on each of them. We need some further definition. Given 
a maximal split component G' of a split pair {s, t} of G* , a vertex w G G' 
properly belongs to G' if w ^ s, t. Given a split pair {s, t} of G* , a maximal 
split component G' of {s, t} is internal if neither u nor v (the poles of G*) 
properly belongs to G', external otherwise. A maximal split pair {s,t} of 
G* is a split pair of G* that is not contained into an internal maximal split 
component of any other split pair {s', t'} of G*. Let {u±, v±}, . . . , {uk, v^} 
be the maximal split pairs of G* (k > 1) and, for i = 1, . . . , k, let Gi be 
the union of all the internal maximal split components of {uj, v{\. Observe 
that each vertex of G* either properly belongs to exactly one Gi or belongs 
to some maximal split pair {u{, Vi}. Node p, is an R-node. Graph skel(//) is 
the graph obtained from G* by replacing each subgraph Gi with the virtual 
edge e, between Ui and Vi. The decomposition recurs on each Gi with \i as 
parent node and with {uj, Vi} as pair of vertices. 



For each node fx of T, the construction of skel(/i) is completed by adding a 
virtual edge (u, v) representing the rest of the graph. 

The SPQR-tree T of a graph G with n vertices and m edges has m Q-nodes 
and 0(n) S-, P-, and R-nodes. Also, the total number of vertices of the skeletons 
stored at the nodes of T is 0(n). Finally, SPQR-trees can be constructed and 
handled efficiently. Namely, given a biconnected planar graph G, the SPQR-tree 
T of G can be computed in linear time [7, 8, 12]. 

2.3 Book Embedding 

A book embedding of a graph G = (V,E) consists of a total ordering of the 
vertices in V and of an assignment of the edges in E to pages, in such a way that 
no two edges (a, b) and (c, d) are assigned to the same page if a -< c -< b -< d. 
A k-page book embedding is a book embedding using k pages. A partitioned 
k-page book embedding is a /c-page book embedding in which the assignment 
of edges to the pages is part of the input. In the special case when k = 2, we 
call the problem PARTITIONED 2-PAGE BOOK EMBEDDING (P2BE). Hence, 
an instance of the PARTITIONED 2-PAGE BOOK EMBEDDING problem is just a 
graph G(V, Ei, E2), whose edges are partitioned into two sets E\ and E2, the 
pages are just two, and the edges of E\ are pre-assigned to page 1 and the edges 
of E2 are pre-assigned to page 2. We say that the edges of E\ (of Ei) are red 
(blue) edges. 

2.4 Eulerian Tour 

Let G be a directed planar embedded graph. A directed cycle of G is a Eulerian 
tour if it traverses each edge exactly once. Consider a vertex v of G and let 
(vi, v), (v, V2), (v, vs), and (t>4, v) be four edges incident to v appearing in this 
order around v in the given embedding. If a Eulerian tour contains egdes (v\ , v), 
(v, Vs), (t>4, v), and (v, v%) in this order then it is self-intersecting. 

3 A Partitioned 2 -page Book-Embedding Testing Algorithm 

In this section we describe an algorithm that, given an instance of P2BE, decides 
whether it is positive and, in case it is, constructs a book embedding of the input 
graph such that each edge is drawn on the page it is assigned to. The algorithm is 
the one proposed in [13]. However, substantial modifications have been applied 
to implement it. Part of them aim at simplifying the algorithm, while others at 
decreasing the value of some constant factors spoiling the efficiency. Further, 
some steps that are described at high level in [13] are here detailed. The main 
differences with [13] are highlighted throughout the paper. 



(a) (b) (c) (d) 

Fig. 1. (a) A disjunctive and splitter-free embedding of a graph, (b) The corresponding green 
graph, (c) An r-rimmed embedding of a graph G. (d) An embedding of G that is not r-rimmed. 



Let G(V, E\,E2) be an instance of problem P2BE. We say that the edges 
of E\ (of E2) are red {blue) edges. As pointed out in [13], the cases in which G 
is disconncted or simply connnected can be easily reduced to the case in which 
G is biconnected, in the sense that G admits a P2BE if and only if all the bicon- 
nected components of G admit a solution. In fact, simply connected components 
can just be placed one after the other on the spine of the book embedding, while 
biconnected components need to be connected through their cut-vertices. How- 
ever, it is easy to see that if a biconnected component admits a book embedding, 
then it admits a book embedding in which the cut-vertex connecting it to its 
parent component in the BC-tree is incident to the outer face. Namely, such a 
book embedding can be obtained by circularly rotating the vertices on the spine. 
Hence, it is always possible to merge the biconnected components on the spine 
through their cut-vertices. Hence, we limit the description to the case in which 
G is biconnected. Moreover, we assume that both E\ and E2 are not empty, 
since a graph with only red (blue) edges is a positive instance if and only if it is 
outerplanar, which is testable in linear time. 

The algorithm is based on a characterization proved in [13] stating that an 
instance admits a solution if and only if G admits a disjunctive and splitter- 
free planar embedding (see Fig. 1(a)). An embedding is disjunctive if for each 
vertex v £ V all the red (blue) edges incident to v appear consecutively around 
v. Notice that, in the upward planarity literature, disjunctive embeddings are 
often called bimodal[ll]. A splitter is a cycle C composed of red (blue) edges 
such that both the open regions of the plane determined by C contain either a 
vertex or a blue (red) edge. An embedding is splitter-free if it has no splitter. 
The first part of the algorithm, that is based on the SPQR-tree decomposition 
of G and whose details are in Sections 4 and 5, concerns the construction of an 
embedding of G satisfying these requirements, if it exists. Otherwise, G does 
not admit any solution. 

Once a disjunctive and splitter-free embedding r of G has been computed, 
an auxiliary graph G* , called green graph, is constructed starting from r. Then, 
as proved in [13], a P2BE of G can be constructed by computing a non-self- 



intersecting Eulerian tour on G* and by placing the vertices of V on the spine 
in the order they appear on such a Eulerian tour. 

Graph G* is a directed graph whose vertices are the vertices of V plus a 
vertex for each face of r. See Fig. 1(b). Edges of G* are determined as follows. 
For each vertex v of G incident to at least one red edge and one blue edge, 
consider each face / incident to v such that v is between a red edge e\ and a 
blue edge e2 on /. If e\ immediately precedes e2 in the clockwise ordering of 
the edges around v, then add to G* an oriented edge (v, /), otherwise add an 
oriented edge (/, v). For each vertex w of G* incident only to red (blue) edges, 
consider a face /' incident to w that contains at least one blue (red) edge. Since 
r is splitter-free, such face exists. Then, add directed edges (w, /') and (/', w). 
Note that, by construction, G* is a bipartite plane digraph, every vertex v of V 
has degree 2 in G* , namely v is incident to exactly one entering and one exiting 
edge, and each vertex / corresponding to a face of r has even degree, namely 
the number of edges entering / equals the number of edges exiting /, and such 
edges alternate around /. From this and from the fact that the underlying graph 
of G* is connected, as pointed out in [13], it follows that G* contains a Eulerian 
tour. 

In the following we show that the alternation of entering and exiting edges 
around each vertex ensures the existence of a non-selfintersecting Eulerian tour, 
as well. In order to do that, we describe an algorithm that, given a disjunctive 
and splitter-free embedding and the corresponding green graph G* , computes a 
non-self-intersecting Eulerian tour of G*. 

Given a plane embedded graph and an outer face /, we call boundary the set 
of (possibly non-simple) cycles composed of edges that are incident to /. We 
proceed on the green graph G* as follows. Starting from any outer face / we 
iteratively remove at each step i the edges of the boundary Bi, thus identifying 
a new outer face and a new boundary, until the graph is empty. On the cycles 
belonging to the extracted boundaries a hierarchical relationship is defined as 
follows. Given two consecutive boundaries and B^+i, a cycle Cj of B^ is 
the father of a cycle Ck of Bh+i if Cj and Cfe share a vertex. This hierarchy 
can be easily represented by a tree, which we call the boundaries tree, whose 
nodes are the cycles of the boundaries and whose root is the cycle representing 
the outer face of G*. Given the alternance of outgoing and incoming edges on 
the nodes of G* , it is easy to see that every connected component of a boundary 
Bi is a directed cycle. A Eulerian non-self-intersecting tour of G* is obtained 
by visiting every cycle of the boundaries according to its orientation in the order 
induced by a DFS visit of the boundaries tree. Namely, starting from an edge of 
the cycle that is the root of the boundary tree, we construct the tour by following 
the orientation of the edges. When a node v of degree greater than 2 is encoun- 



tered coming from an oriented edge (u, v) of a cycle, we start visiting its child 
cycle by following the edge (v, w) following (u, v) in the clockwise order of 
the edges around v. Note that, because of the alternance of entering and exiting 
edges, edge (v, w) is directed from v to w. The same happens when the visit of 
the child is finished and the visit of the father continues. Hence, intersections in 
the Eulerian tour are always avoided. 

From the above discussion it follows the claimed statement that the de- 
scribed algorithm computes a P2BE of (V, E\,E2), if any such a P2BE exists. 

4 Computing a Disjunctive and Splitter-Free Embedding 

Let G(V, Ei, E2) be a biconnected planar graph. We describe an algorithm to 
compute a disjunctive and splitter-free embedding of G, if any such an embed- 
ding exists, consisting of two preprocessing traversals of the SPQR-tree T of G 
and of a final bottom-up traversal to compute the required embedding. 

Let /i be a node of T ■ According to [13], a virtual edge e of skel(/i) is an r- 
edge (a b-edge) if there exists a path in pert(/x) between the poles of p, composed 
of red edges (of blue edges). If e is both an r-edge and a b-edge, it is a br-edge. 

Consider a cycle C = e\, . . . , e q in skel(/i) composed of edges of the same 
color, say r-edges. If C is a splitter in every embedding of skel(/x), then a splitter 
is unavoidable. However, even if there exists an embedding of skel(/x) such that 
C is not a splitter, then a cycle in pert(//) passing through the pertinent graphs 
of ei, . . . , e q could still be a splitter (since e±, . . . , e q are r-edges, there exists 
at least one red cycle C in pert(/i)). Consider any node v corresponding to 
a virtual edge ej and the path p v {C) between the poles of v that is part of C . 
Intuitively, in order for C not to be a splitter, we should construct an embedding 
of pert(z/) in which p u {C) is on the outer face. Actually, not all the vertices of 
p u (C) have to be on the outer face, since red chords might exist in p u {C) 
(that is, red edges connecting vertices not consecutive in p v {C')), separating 
some vertex of p v (C) from the outer face, as in this case such chords would be 
internal to C, and this does not make it a splitter. On the other hand, if p v {C) 
has a blue edge or a vertex (even if this vertex belongs to another path between 
the poles composed of red edges) on both its sides, then C becomes a splitter. 
In analogy with [13], where the same concept was described with a slightly 
different definition, we say that an embedding of pert(!/) in which each path 
between the poles composed of red edges (of blue edges) has only red edges 
(blue edges) on one of its sides is r-rimmed (is b-rimmed). Figs. 1(c) and (d) 
show an r-rimmed and a non-r-rimmed embedding, respectively. Note that an 
embedding could be at the same time both r- and b-rimmed, with the red and the 
blue paths on different sides of the outer face. 
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Fig. 2. Parallel virtual edges are sketched with rectangles colored according to their poles, (a) An 
r-rimmed embedding forces an RBR color-pattern on a pole, (b) A color-pattern BR or RB on a 
pole forces either an RBR or a BRB on the other pole, (c) An R-node. Virtual edges representing 
Q-nodes are thin, (d) The corresponding auxiliary graph Oi. (e) A splitter that is not a rigid- 
splitter, (f) Disjunctiveness constraints on nodes ei and ei determine a splitter (ei, e2, e^, e^). 

The existence of an r-rimmed (b-rimmed) embedding is necessary only for 
each node fi such that there exists a cycle C of red (blue) edges traversing both 
fx and its parent. However, the existence of C is not known when processing [i 
during a bottom-up visit of T. Thus, we perform a preprocessing phase to decide 
for each node \i whether any such cycle C exists. In this case, /x is r-joined (b- 
joined). Hence, when processing fj,, we know whether it is r-joined (b-joined) 
and, in case, we inductively compute an r-rimmed (b-rimmed) embedding. 

Concerning disjunctiveness, for each vertex w of skel(/i) we have to check 
whether the ordering of the edges around w determined by the embedded per- 
tinent graphs of the child nodes incident to w makes it disjunctive. In order to 
classify the possible orderings of edges around the poles of a node we define, in 
analogy with [13], the color-pattern of a node /iona vertex v as the sequence 
of colors of the edges of pert(/i) incident to v. Namely, the color-pattern of fi 
on v is one of R, B, RB, BR, RBR, BRB. Note that, if the color-pattern is 
either R or B, then it is the same in any embedding. Otherwise, it depends on 
the chosen embedding. Hence, it might be influenced by the fact that the em- 
bedding needs to be r- or b-rimmed (see Fig 2(a)) and by the need of a particular 
color-pattern on the other pole (see Fig 2(b)). Thus, a color-pattern either RBR 
or BRB could be forced on a pole u of fi although an RB or a BR pattern 
would be possible as well. Another factor influencing the color-pattern on u is 
the presence of red or blue edges incident to u in the pertinent of the parent v of 
[i. In fact, if u has color-pattern RBR (BRB) and there is a blue (red) edge in 
pert(i/) incident to u, then u is not disjunctive. Thus, in the preprocessing phase 
we also determine two flags for each pole u of \x, stating whether v contains 
at least one red (blue) edge incident to u. Hence, when processing fi, we know 
whether it is admissible to have an RBR (a BRB) color-pattern on its poles. 

Hence, after the preprocessing phase, we can assume to know for each node 
H the following information: 

1 . two flags stating whether \i is r-joined and whether it is b-joined; 




2. for each pole u of /j,, two flags stating whether the parent v of li contains 
at least one red edge and whether it contains at least one blue edge, respec- 
tively, incident to u. 

The two information obtained in the preprocessing can be properly com- 
bined when processing a node to decide whether an embedding satisfying all 
the constraints exists, as described is Section 5. If it is not the case, we state that 
the instance is negative, while in the case that at least one of such embeddings 
exists, we can arbitrarily choose one of them, without the need of carrying on a 
multiplicity of embeddings. This is one of the most crucial differences between 
our implementation and [13]. In fact, even if they perform a preprocessing 
to determine whether a node is r-joined (b-joined), they do not exploit it for 
disjunctiveness, and have to consider at each step all the possible embeddings 
determining different color-patterns on the two poles. Of course, as the number 
of color-patterns is bounded by a constant, this does not affect the asymptotic 
complexity, but our solution noticeably improves on the execution times. Also, 
they deal with constraints given by the r-joinedness (b-joinedness) and by the 
disjunctiveness in two different steps. In our case, instead, instances that are 
negative due to disjunctiveness are recognized much earlier. 

The preprocessing consists of a bottom-up and a top-down traversal of T, 
that we describe in the following. The bottom-up traversal computes some in- 
formation on each node, which are then aggregated in the top-down traversal to 
efficiently compute the needed information on the parent of each node. 

In order to determine which are the r- and the b-joined nodes, in the bottom- 
up traversal we compute for each node whether its skeleton (excluding the vir- 
tual edge representing the parent) contains a path between its poles composed of 
r-edges (b-edges). Then, in the top-down traversal we transmit this information 
from each node to its children, namely all and only the children that are part 
of a cycle composed of r-edges (b-edges) in the skeleton of a node are r-joined 
(b-joined). 

In order to determine which are the nodes whose parent has at least a red (a 
blue) edge incident to a pole u, we determine for each node ll in the bottom-up 
traversal whether it contains a red (blue) edge incident to u, and in case it does, 
we add 1 to a counter associated with u and the parent of li. Then, during the 
top-down traversal we inductively compute the information on each node li, we 
accordingly update the counter associated with u and ll for each child node of li, 
and we state that the parent of a child node v of ll has a red (blue) edge incident 
to u either if the value of the counter is at least 2 or if it is 1 and v has no red 
(blue) edge incident to u. 

In the next section we describe the final bottom-up traversal of T which 
computes a disjunctive and splitter-free embedding of G, if it exists. 



5 SPQR-tree Algorithm 



When considering a node fiofT with children ui, . . . , v^, exploiting the infor- 
mation resulting from the preprocessing and the information inductively com- 
puted for fx, . . . , Vk, we check whether /i admits a splitter-free and disjunc- 
tive embedding and compute the following: (i) if \i is r-joined (b-joined), an 
r-rimmed (a b-rimmed) embedding; and (ii) the color-patterns of the poles of p. 

In the base case, \i is a Q-node. Suppose that skel(/i) is an r-edge, the other 
case being analogous. If fi is r-joined, every embedding of skel(/i) is r-rimmed. 
Further, the color-pattern on the poles is R in any embedding of skel(yu). 

Suppose that fi is an R-node. Since skel(/i) is triconnected, it has one planar 
embedding, up to a flip. Hence, if there is a splitter in skel(^), then it is unavoid- 
able. Hong and Nagamochi call such splitters rigid-splitters. In order to test the 
existence of such splitters, for each set Ei, i = 1, 2, we construct an auxiliary 
graph Oi starting from skel(/i). See Figs. 2(c) and (d). We describe the construc- 
tion for Ei, the other case being analogous. Initialize 0± =skel(/i). Subdivide 
each virtual edge of skel(^) (including the one representing the parent) with a 
dummy vertex, except for the r-edges corresponding to Q-nodes. Then, for each 
dummy vertex subdividing a virtual edge that is not an r-edge, remove one of 
its incident edges without modifying the embedding. Finally, check whether the 
obtained embedding of 0\ is an outerplane embedding, that is, all the vertices 
of 0\ are on the same face. This check is performed by iterating on all the faces 
of the embedded graph 0\ and by checking whether there exists one containing 
all the vertices. Note that this step can be performed in linear time, since each 
vertex of degree d is examined at most d times and since the sum of the degrees 
of the vertices of a graph is twice the number of edges, which is 0(n). In [13] 
this step is performed by constructing a variant of the green graph and checking 
whether it is connected. Even if the time complexity of the two approaches is 
basically the same, we find that our approach is easier to implement and slightly 
more efficient, since 0\ does not need to be constructed, but can be obtained by 
flagging the edges of skel(^/). 

Note that, for each cycle composed of r-edges (b-edges) in skel(/i) that is 
not a rigid-splitter, all the nodes composing it inductively admit an r-rimmed 
(b-rimmed) embedding. Hence, it suffices to flip them in such a way that their 
red (blue) border is turned towards the red (the blue) outerplanar face. However, 
if each of them has an embedding that is both r-rimmed and b-rimmed, the red 
and the blue outerplanar faces coincide and it is not possible to flip the nodes 
properly, which implies that a splitter exists in the embedding. See Fig. 2(e). 
This type of splitter seems to have gone unnoticed in [13], where flips imposed 
by cycles of r- and b-edges are considered independently. 



We deal with disjunctiveness constraints. We observe some straightforward 
properties of the color-patterns of the nodes incident to the same vertex w of 
skel(/i). (i) At most two nodes have color-pattern different from R and B. (ii) If 
one node has color-pattern RBR (BRB), then all the other nodes have color- 
pattern R (B). Hence, since each vertex has degree at least 3 in skel(/i), at least 
one node v incident to w exists with color-pattern either R or B. Thus, starting 
from v, we consider all the nodes incident to w in clockwise order and greedily 
decide a flip based on the current color. If more than two changes of color are 
performed, then G does not admit any disjunctive embedding. If exactly one 
node v has color-pattern different from R or B and all the other nodes have 
color-pattern R (B), then the flip of v is not decided at this step. Also, the flip 
is not decided for the nodes having color-pattern R or B. 

Disjunctiveness and splitter-free constraints might be in contrast. See Fig. 2(f). 
We can efficiently determine such contrasts by flagging the nodes that need to be 
flipped and, in case such contrasts exist, state that the instance is negative. This 
check is not described in [13], where possible contrasts between disjunctive and 
splitter-free constraints are noticed for P-nodes but not for R-nodes. 

The color-patterns of the poles and, if needed, an r-rimmed (a b-rimmed) 
embedding of pert(/i) are computed by considering the information on the parent 
node, the color-patterns of the virtual edges incident to the poles, and the r- 
rimmed (b-rimmed) embedding of the children. 

Suppose that \i is an S-node. Since skel(/i) is a cycle containing all the 
virtual edges, even if such a cycle is composed of edges of the same color, then 
it is not a splitter. Namely, even if there exist both a red and a blue cycle passing 
through all the children of [i, such nodes can be flipped so that the red and the 
blue borders are turned towards the two faces of skel(^). 

Concerning disjunctiveness constraints, if two children both incide on a ver- 
tex u of skel(/i) with color-pattern either BR or RB, then they have to be flipped 
in such a way that the red edges (and hence the blue edges) are consecutive 
around u. In all the other cases, the relative flip of the two children incident to 
u is not fixed by their color-patterns. If there exists at least a vertex u with this 
property, we say that /i admits two different semi-flips. Intuitively, this means 
that the color-pattern of a pole is independent of the one on the other pole, since 
they depend on flips performed on two different subsets of children of \i. 

Note that in an S-node no contrast between splitter-free and disjunctiveness 
constraints are possible, since flipping the r-rimmed embeddings towards the 
same face implies placing the red edges consecutive around u. Hence, no nega- 
tive answer can be given during the processing of an S-node. 



The color-pattern on each pole is the color-pattern of the unique node inci- 
dent to it, while an r-rimmed (b-rimmed) embedding is obtained by concatenat- 
ing the r-rimmed (b-rimmed) embeddings of the children. 

Suppose that fi is a P-node. In order for a splitter-free embedding to exist, 
the following must hold: (i) There exist at most 3 r-edges (b-edges); if they are 
3 then one is a Q-node. (ii) There exist at most 2 virtual edges that are both 
r-edges and b-edges; if they are 2 then there exists only another virtual edge and 
it is a Q-node. When such conditions do not hold, the r-edges (b-edges) induce 
a splitter in every embedding of the P-node. 

On the other hand, in order for a disjunctive embedding to exist, the follow- 
ing must hold: (i) if there exists a virtual edge with RBR (BRB) color-pattern 
on a pole, then all the other edges have color-pattern R (B) on that pole; (ii) there 
exist at most two virtual edges with color-pattern RB or BR on a pole. When 
these conditions do not hold for a pole u, in every embedding of G there exist 
more than two color changes in the clockwise ordering of edges incident to u, 
that is, there exists no embedding that makes u disjunctive. 

Consider a child node v\ having color-pattern either R or B on both poles, 
say R on pole u and B on pole v, and consider another child node i>2 having 
color-pattern Ronu and B on v. Nodes ui and V2 can be considered as a single 
node v* with color-patterns R and B on the two poles. When the permutation of 
the P-node has been computed, v* is replaced by v\ and V2 . This operation re- 
duces the number of virtual edges to at most 8, namely at most 4 groups of nodes 
having either R or B on both poles plus at most 2 nodes with color-pattern dif- 
ferent from R and Bona pole and at most 2 nodes with color-pattern different 
from R and B on the other pole. Note that the parent cannot be grouped, since its 
color-patterns are unknown at this stage. In [13] this fact is exploited to search 
an embedding with the desired properties by exhaustively checking all permu- 
tations, i.e., with a brute-force approach. However, even if the time complexity 
is aympthotically linear, this yields a huge number of cases, namely 8! * 2 8 
combinations, i.e., all permutations of 8 edges multiplied by all flip choices. 

Hence, our implementation uses a different approach in order to search 
into a much smaller space. Namely, consider any color-pattern, say RBR, and 
map it to a linear segment of fixed length, partitioned into three parts R, B, 
R, by two points that represent the two changes of color R — B and B — R. 
Such points are identified by a unidimensional coordinate p along the segment. 
Given two color-patterns, their segments, and a separating point for each of 
them, with coordinates pi and p2, respectively, any of the following condi- 
tions can hold: (i) p\ < p2, (ii) pi = P2\ (hi) Pi > P2- See Figure 3(a). We 
call alignment of two color-patterns each combinatorial possibility obtained by 
exhaustively making conditions (i)-(iii) hold for all pairs of separating points 
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Fig. 3. (a) All possible alignments of a pair of RB color-patterns, (b) Correspondence between 
an alignment of a pair of color-patterns and a sequence of virtual edges of a P-node. 

of their segments. An alignment of two color-patterns P\,P-z uniquely corre- 
sponds to a sequence of virtual edges whose color-patterns compose Pi and 
P2 on the two poles. Such sequence makes both poles disjunctive by construc- 
tion. See Fig. 3(b). Our approach exploits this fact by exhaustively enumerat- 
ing all alignments of all pairs of color-patterns. The result is the set L con- 
taining all and only the disjunctive edge permutations of a generic P-node. 
L contains exactly 180 elements. The pseudocode of the algorithm Gener- 
ate_Admissible_Solutions_for_P-nodes performing the enumeration that 
constructs L is given in Algorithm 1 . 

Since the virtual edges of a P-node have a disjunctive permutation if and 
only if they can be disposed in the same sequence as an element in L, a disjunc- 
tive embedding can be found, if it exists, by a brute force search across the 180 
elements of L, an impressive improvement with respect to the algorithm in [13]. 

As the parent node could not be grouped with other nodes, it could impose 
some additional constraints on the permutation to find that forbid permutations 
having color-patterns RBR or BRB and that require any r-rimmed (b-rimmed) 
node to be either the first or the last, if the P-node is r-joined (b-joined). 

The whole P-node algorithm must be repeated for every possible choice of 
semi-flip for the virtual edges admitting it. However, at most two such virtual 
edges can exist, since they have color-patterns RB or BR on both poles. Hence, 
the algorithm must be repeated up to 4 times. 

6 Experimental Results 

In this section we describe the experimental tests performed to check correct- 
ness and efficiency of our implementation. When performing experiments a cru- 
cial aspect is to have at disposal a representative set of negative and positive 
instances. Negative instances have the main role of checking the correctness, 
while positive instances are both used to check the correctness and to test the 
performance in a complete execution, without being influenced by early recog- 
nition of negative instances. We constructed the former set using ad-hoc exam- 
ples, conceived to stress all the steps of the algorithm. On the other hand, in or- 



Algorithm 1 Generate_Admissible_Solutions_for_P-nodes 



1. L <— empty list of permutations of virtual edges 

2. S <- list of color-patterns: R, B, RB, BR, RBR, BRB 

3. for all (Ti€5do 

4. for all (72 G S do 



5. for all alignment of ai and 02 do 

6. Z <— list of pairs of colors, where each pair is composed of a color of o\ and of a 
color of 02. Elements of the list are determined by discretizing the alignment. Note 
that each alignment determines at most nine list elements. 

7. P <— empty list of virtual edges 

8. for all z G Z do 

9. append to P a new edge p with color-patterns on its poles G {R, B} corresponding 
to the colors of z 

10. if the color-patterns of the poles of p are either both R or both B then 

1 1 . make p r-rimmed or b-rimmed depending on whether the color-patterns are R 
01B 

12. end if 

13. if p is the first or the last element of Z then 

14. flip p in such a way that the r-rimmed (b-rimmed) path is towards outside 

15. end if 

16. end for 

17. for all p G P do 

18. if the first color of the color-pattern of p on a pole is different from the last color 
of the color-pattern of the edge preceding p in P on the same pole then 

19. insert a new edge p' preceding p with color-pattern RB or BR on the considered 
pole 

20. end if 

21. if the color-patterns of p' on the two poles have the same first (last) color then 

22. make p' r-rimmed or b-rimmed 

23. end if 

24. end for 

25. D <— edges that have color-pattern R or B on both poles and multiple instances in P 

26. if size(D) = then 

27. append(P, L) 

28. else 

29. for all p G D do 

30. for all instance pi of p in P do 

31. P' «— copy of P with only instance pi of p 

32. append(P',L) 

33. end for 

34. end for 

35. end if 

36. end for 



37. end for 

38. end for 

39. return L 



der to obtain a suitable set of positive instances, we used random generation. Un- 
fortunately, to the best of our knowledge, no graph generator is available to uni- 
formly create graphs with a P2BE. Hence, we devised and implemented a graph 
generator, whose inputs are a number n of vertices and a number m < 3n — 6 of 
edges. The output is a positive instance of P2BE selected uniformly at random 
among the positive instances with n vertices and m edges. 

The generator works as follows. First, we place n vertices v\, . . . , v n on the 
spine in this order. Then we insert, above (below) the spine, red (blue) dummy 
edges (vijVz), ■ ■ ■ , (v n -i,v n ), and (vi,v n ). In this way we initialize the two 
pages with two faces (v±, . . . , v n ) composed of red and of blue dummy edges, 
respectively. Observe that we inserted multiple dummy edges. Dummy edges 
will be removed at the end. Second, we randomly select a face / with at least 
three vertices, selected with a probability proportional to the number of candi- 
date edges that can be added to it. Then, an edge (u, v) is chosen uniformly at 
random among the potential candidate edges of /. Edge (u, v) is added to / by 
either splitting / or substituting a dummy edge of / with a "real" edge. Edge 
(u, v) is colored red or blue according to the color of the edges of /. If (u, v) is 
red (blue) we check if there exists a blue (red) face that contains both u and v 
and remove (u, v) from the candidate edges of that face. We iteratively perform 
the second step until m is reached and at the end we remove the dummy edges 
that have not been substituted by a "real" edge. Observe that in this way we do 
not generate multiple edges and that the generated graphs are not necessarily 
connected. 

We generated three test suites, Suite 1, 2, and 3, with m = 2n, m = 2.5n, 
and m = 3n — 6, respectively. For each Suite, we constructed ten buckets of 
instances, ranging from n = 10, 000 to n = 100, 000 with an increment of 
10, 000 from one bucket to the other. For each bucket we constructed five in- 
stances with the same parameters n and m. The choice of diversifying the edge 
density is motivated by the wish of testing the performance of the algorithm on 
a wide variety of SPQR-trees, with Suite 3 being a limit case. 

The algorithm was implemented in C++ with GDToolkit [6]. The OGDF 
library [4] was used to construct the SPQR-trees. We used GDToolkit because 
of its versatile and easy-to-use data structures and OGDF to construct SPQR- 
trees in linear time. 

Among the technical issues, the P-node case required the analysis of a set 
of cases that is so large to create correctness problems to any, even skilled, 
programmer. Hence, we devised a code generator that, starting from a formal 
specification of the constraints, wrote automatically the required C++ code. For 
performing our experiments, we used an environment with the following fea- 
tures: (i) CPU Intel Dual Xeon X5355 Quad Core (since the algorithm is se- 
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Fig. 4. Execution times of the generator for generating the three suites. The a>axis represents the 
number of vertices of the bucket, while the y-axis represents the average execution time on the 
instances in the bucket. 
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Fig. 5. The a;-axis represents the number of vertices of the bucket, while the y-axis represents the 
average number of components of each bucket. 



quential we used just one Core) 2.66GHz 2x4MB 1333MHz FSB. (ii) RAM 
16GB 667MHz. (iii) Gentoo GNU/Linux (2.6.23). (iv) g++ 4.4.5. 

Figs. 4(a)-(c) show the execution times of the generator for generating the 
three suites. 

Before giving the execution times of the algorithm on the generated in- 
stances, we show some charts describing the structure of such instances, both 
in terms of connectivity and in terms of the complexity of the corresponding 
SPQR-trees. 

Figs. 5(a)-(c) show the number of connected (including isolated vertices) 
and biconnected (including single edges) components in the test suites. 
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Fig. 6. Number of SPQR-tree nodes in the three test suites. The y-axis represents the average 
number of SPQR-tree nodes of each bucket. 
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Fig. 7. Total execution time of the algorithm on the three test suites. The y-axis represents the 
average execution time of the algorithm on the instances in the bucket. 



Figs. 6(a)-(c) show the number of SPQR-tree nodes in the three test suites. 
Note that the large amount of P-nodes in Suites 1 and 2 puts in evidence how 
crucial has been in the implementation to optimize the P-nodes processing. 

Then, we give the execution times of the algorithm on such instances and an 
analysis of them from several points of view. 

Fig. 7(a)-(c) show the total execution times for the three suites. These mea- 
surements include the time necessary to decompose the graphs in their con- 
nected, biconnected, and triconnected components. The algorithm clearly shows 
lineal - running times, with very little differences among the three suites. 

Figs. 8(a)-(c) show the execution times of the main algorithmic steps for the 
three suites, namely (i) the total time spent to process biconnected components, 
(ii) the time spent to deal with the SPQR-trees (excluding the time to create 
them), and (iii) the time spent to create the green graphs and to find the Eulerian 
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Fig. 8. Execution times of the algorithm for the biconnected components of the three test suites. 
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Fig. 9. The execution times of the four algorithmic substeps of the step that deals with the SPQR- 
trees (excluding creation) for the three suites. 



tours. Beside remarking the linear running time, these charts show how the time 
spent on biconnected components is distributed among the two main algorithmic 
steps. 

Figs. 9(a)-(c) show the execution times of the four algorithmic substeps of 
the step that deals with the SPQR-trees (excluding creation) for the three suites. 
Namely, the five curves show: (i) the time to deal with the SPQR-trees (ex- 
cluding the time spent to create them), (ii) preprocessing bottom-up phase, (iii) 
preprocessing top-down phase, (iv) the bottom-up skeleton embedding phase, 
and (v) the pertinent graph embedding phase. 

Figs. 10(a)-(c) show the execution times of the four algorithmic substeps 
of the step that deals with the green graph for the three suites. Namely, the five 
curves show: (i) the time to deal with the green graph, (ii) creation phase, (iii) 
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Fig. 10. The execution times of the four algorithmic substeps of the step that deals with the green 
graph for the three suites. 
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Fig. 11. The time spent to deal with the different types of SPQR-tree nodes (excluding creation) 
for the three suites. 



decomposition phase, (iv) Eulerian tour finding phase, and (v) Eulerian tour 
visiting phase. 

Figs. 1 l(a)-(c) show the time spent to deal with the different types of SPQR- 
tree nodes (excluding creation) for the three suites. Namely, the four curves 
show: (i) the time for S-nodes, (ii) the time for P-nodes, (iii) the time for Q- 
nodes, and (iv) the time for R-nodes. 



7 Conclusions 

We described an implementation of a constrained version of the 2-page book 
embedding problem in which the edges are assigned to the two pages and the 
goal is to find an ordering of the vertices on the spine that generates no crossing 



on each page. The implemented linear time algorithm is the one given in [13], 
with several variations aimed at simplifying it and at improving its performance. 

We performed a large set of experimental tests on randomly generated in- 
stances. From these experiments we conclude that the original algorithm, to- 
gether with our variations, correctly solves the given problem, and that its per- 
formance are pretty good on graphs of medium-large size. 
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